Interaction and excitonic insulating transition in graphene 
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The strong long-range Coulomb interaction between massless Dirac fermions in graphene can 
drive a semimetal-insulator transition. We show that this transition is strongly suppressed when 
the Coulomb interaction is screened by such effects as thermal fluctuation, doping, disorder, and 
finite volume. It is completely suppressed once the screening factor /i is beyond a threshold 
even for infinitely strong coupling. However, such transition is still possible if there is an additional 
strong contact four-fermion interaction. The differences between screened and contact interactions 
are also discussed. 
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The low-energy elementary excitations in undoped 
graphene are massless Dirac fermions. Their spectral 
and transport properties are quite unusual and have at- 
tracted intense investigations in the past several ye&rs^. 
For a clean undoped graphene, the density of states 
(DOS) N{ijj) vanishes linearly near the Dirac point. 
As a result, the Coulomb interaction between massless 
Dirac fermions is essentially unscreened, in sharp con- 
trast to the electron system with parabolic dispersion. 
The unscreened, long-range Coulomb interaction was 
shown to be responsible for many anomalous behaviors 
of graphcneM^^^i^^OOdi . 

At the strong coupling regime, the long-range Coulomb 
interaction can open a finite mass gap for the Dirac 
fermion, which then drives a phase transition from the 
semimetal state to an insulator state. This transition is 
realized by forming stable particlc-holc pairs and usu- 
ally named as excitonic semimetal-insulator (SM-IN) 
transition^. Recently, this kind of phase transition has 
been studied by nonperturbative Dyson-Schwinger (DS) 
equation approacb^i^, renormalization group method^, 
and lattice simulationai^ii^. The SM-IN transition was 
found in graphene for strong Coulomb coupling and small 
fermion flavori^. The effects of finite temperature and 
external magnetic field were also considered^. 

Although being of remarkable interests, the predicted 
SM-IN transition (in zero magnetic field) has not yet 
been unambiguously observed in experiments. In this 
paper, we discuss the effects that can potentially prevent 
the appearance of this SM-IN transition. First of all, 
it should be emphasized that such transition can take 
place only for strong, poorly screened Coulomb inter- 
action. Generically, there are two critical parameters: 
critical dimensionless coupling strength Ac and critical 
fermion flavor Nc- SM-SI transition is possible only when 
N < Nc and A > Ac. Once the long-range Coulomb in- 
teraction is screened by some physical effects, there will 
be an effective screening factor /i, which is expected to 
increase Ac and reduce Nc. This can be understood by 
noting the important fact that SM-SI transition realized 
by forming fermion-antifermion pairs is a genuine low- 
energy phenomenon. From the experience in QED3, the 
long-range nature of gauge interaction plays the crucial 



role in generating the dynamical mass gap for initially 
massless Dirac fermiona^^. A finite gauge boson mass 
rapidly reduces the critical fermion flavor to below the 
physical value 3ii. In the present case, there is a similar 
suppressing effect once the long-range Coulomb interac- 
tion is screened for some reason. The opening of excitonic 
gap requires that the Coulomb interaction is sufficiently 
strong at low-momentum region. However, the screening 
factor fi suppresses the contribution from small momenta 
significantly. Obviously, this kind of pairing instability is 
markedly different from the conventional BCS-type pair- 
ing formation, which is caused by arbitrary weak attrac- 
tive force between electrons. 

In realistic graphene samples, the critical behavior of 
SM-IN transition can be influenced by the following rea- 
sons; thermal fluctuation; doping; disorder; finite sample 
volume. Each of them can generate an effective screening 
factor /i, which could be regarded as an effective photon 
mass. We study their effects on critical strength Ac and 
critical flavor Nc by solving the corresponding gap equa- 
tion, and show that a growing /j, significantly increases 
Ac and reduces Nc, both at zero and finite temperatures. 
When n is beyond some threshold /Xc, the excitonic tran- 
sition is completely prohibited, leaving semi-metal as the 
stable ground state. Frequently, some of these effects 
coexist, leading to further suppression of excitonic tran- 
sition. However, even when /i > /ic, we found that the 
excitonic transition can still take place if there is an addi- 
tional strong contact quartic interaction. We also briefly 
discuss the interesting differences between the screened 
Coulomb and contact quartic interactions. 

The total Hamiltonian of massless Dirac fermion H = 
Hq + He is given by 
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Here, we adopt four-component spinor field ip to describe 
the massless Dirac fermion since there is no chiral symme- 
try in the two-component representation. The conjugate 
spinor field is defined as ■0 = V'Wo- The 4x4 7-matrices 
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satisfy the standard Clifford algebra. Although the phys- 
ical fermion flavor is actually iV = 2, in the following we 
consider a large N in order to perform 1/A^ expansion. 
The total Hamiltonian preserves a continuous U(2N) chi- 



ral symmetry 



■0, which will be dynamically 



broken if a nonzero fermion mass gap is generated. 

The free propagator of masslcss Dirac fermion is 
G'o(fcojk) = (7oA:o — WF7'k)^^. The Coulomb interaction 
modifies it to the complete propagator 



G(fco,k) 
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where m(k) denotes the dynamical fermion mass and Ai^2 
the wave function renormalization functions. To the lead- 
ing order in 1/N expansion, the DS integral equation is 
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where the vertex function has already been approximated 
by the bare matrix 79. The nontrivial solution m{p) of 
this equation signals the opening of an excitonic gap. 

In the DS gap equation, V{q) is the Coulomb interac- 
tion function. The bare, unscreened Coulomb interaction 
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has the form Va{q) = in the momentum space. For 
an interacting electron gas, the collective density fluc- 
tuations screen the bare Coulomb interaction Vo(g) to 
V~^((l) = VQ^{q) — ir{q). For ordinary non-rclativistic 
electron gas, the static polarization function Tr{qQ = 0) 
is just the zero-energy DOS, iV(0), which is known to be 
finite. The parameter iV(0) defines the inverse Thomas- 
Fermi screening length. The case for undoped clean 
graphcne is quite different because of the linear disper- 
sion of Dirac fermions. The leading contribution to po- 
larization function is given by 7ro(q) = — ; , , 

It vanishes linearly as q ^ in the static limit q^ = Q, 
so the long-range Coulomb interaction is unscreened. 

Under the approximations described above, the gap 
equation can be written as 
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with interaction function 



V{q) = 
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Here, A1.2 = 1 is assumed and the rescaling Ui?k k, 
Wi?A A is made (such rescaling will be made through- 
out the whole paper). The present problem contains 
two parameters: fermion flavor N and dimensionless 
Coulomb coupling defined as A = g^N /IQvp, where 
gc = e^/eo. The ultraviolet cutoff A is taken to be of 
order lOeV which is determined by ^ with lattice 
constant a ~ 2.46A. Not that no instantaneous approxi- 
mation for the polarization function is made at present. 
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FIG. 1: (a) Relationship between Nc and A at zero tempera- 
ture; (b) Relationship between Nc and A at different temper- 
atures T. Botli are for unscreened Coulomb interaction. 



We solve the nonlinear gap equation using bifurcation 
theory and parameter embedding mctho d^^'^^ for a num- 
ber of fixed values of A. The fermion flavor TV serves as 
the embedded parameter in seeking the bifurcation point. 
The results are shown in Fig.[T](a). It is easy to see that 
the critical flavor TVc is an increasing function of A. For 
A ^ oo, TVc « 3.52; for A = 2, A^^ « 2. 

The above results are valid only for the unscreened 
Coulomb interaction at zero temperature. In realistic 
systems, the long-range interaction could be screened by 
several physical effects, such as thermal fluctuation, dop- 
ing, disorder, and finite volume. If the polarization func- 
tion T:{qQ,q) takes a finite value due to some mechanism 
in the qo = and q ^ limit, then the long-range 
Coulomb interaction becomes short-ranged and 7r(0,0) 
defines the screening factor. Before computing 7r((7o,q) 
by taking each screening effect into account, we now phe- 
nomenologically introduce a single parameter ^ (in unit 
of eV) to model the screened interaction function 



Viq) 
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The advantage of this parameter is that it explicitly mea- 
sures the suppressing effect on the critical behavior due 
to all possible screening mechanisms. If we regard this 
function as the effective interaction strength, then the 
influence of n becomes clear: it eliminates the contribu- 
tion of small momenta to the gap equation Eq. ([3]). But 
remember that the excitonic gap generation is primarily 
determined by the contribution from this region, so it is 
expected that a large /i will destroy SM-IN transition. 
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After solving the gap equation, we found that a growing 
/i leads to increase of critical strength Ac and to decrease 
of critical flavor (see Fig. EJa))- Beyond some critical 
value fj-c, the SM-IN transition is completely prevented, 
even when the dimensionless strength A — > oo. 

The possible screening mechanisms will be discussed in 
order. First of all, the thermal fluctuation will surely re- 
store the chiral symmetry even it is broken by the ground 
state. At finite temperatures, the Matsubara fcrmion 
propagator is 
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where ujn = (2n + l)i:T is the fermion frequency. Here, 
in order to carry out the frequency summation appear- 
ing in the gap equation, we utilize the instantaneous 
approximation^. Under this approximation, the polar- 
ization function can be approximated^^ by 



«^q + crcxp(-^)), (7) 
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with constant c = 161n2/7r. At the limit q ^ 0, the po- 
larization is ~ T, corresponding to the thermal screening 
factor /i. Since other screening effects can coexist with 
thermal fluctuations at finite temperatures, we still in- 
troduce the parameter fj. and write the gap equation as 



-(P,T) = 1/ 



i(k,r) 



87r2 y^k2-hm2(k, T) 
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with the interaction function 
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The results at finite temperatures are rather complex 
since now we have four parameters, N , A, T, /i, each of 
which has a critical value. Their relationships are shown 
in Fig. [ijb) without screening effects = 0) and in 
Fig. [2lb) with screening (the temperature is in unit of 
eV). In Fig. EJb), the Coulomb coupling parameter is 
fixed at A — > oo, and the results for other values of A 
are not shown since they are qualitatively similar. The 
results tell us that the thermal suppression is more im- 
portant than screening effect when /.t has small values 
(< 10~^), but the screening effect eventually becomes 
much more important than thermal effect for larger val- 
ues of /.t. 

The second potential mechanism that can prevent gap 
generation is doping. The Coulomb interaction between 
Dirac fermions is unscreened only when the graphene is 
undoped. When the graphene is shghtly doped, the finite 
carrier density then serves as an effective screening factor 
fi. Then the excitonic gap is expected to open only at or 
very close to the Dirac point. The critical carrier density 



FIG. 2: (a) Dependence of Nc on fj, for different A at zero 
temperature; (b) Dependence of Nc on ^ for different T at 
A — > oo. 



has been discussed previously in^. Recently, the same 
screening effect was emphasized in the study of exciton 
condensate in bilayer grapheneii. At finite chemical po- 
tential /io, the fermion propagator becomes 



C?(iw„,k,^o) 
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Using this propagator, the polarization function can be 
calculated with the result 
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in the zero frequency limit. As q ^ 0, 7r(0, 0, ^o) = 
^/io, which defines the screening factor /i. After per- 
forming the frequency summation, the gap equation has 
the form 
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m(p,r) = — / — -- 
' N J 87r2 
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with the interaction function being 
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Note the chemical potential /ip appears in two places: the 
occupation number and the polarization function. To 
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see the dominant effect of /xq, we first solved the full 
gap equation and show the results in Fig. [3]^ a) (also at 
X —I- oo for comparison). The dependence of Nc on T and 
/io qualitatively resembles that in Fig. [^l^b), but visibly 
exhibits different quantitative behavior: the suppressing 
effect from doping is more prominent at low T than at 
higher T. Despite the details, a large doping makes the 
excitonic insulating state impossible. Then we solved the 
gap equation by ignoring the /Uo-dependence of occupa- 
tion number and found that the results are nearly the 
same as Fig. Ela) (therefore not shown). It seems that 
the screening effect induced by doping plays the domi- 
nant role in suppressing the gap generation. 

Next, we consider the influence of disorders, which 
are unavoidable in graphene samples. The disorders can 
be crudely classified as random mass, random chem- 
ical potential, and random vector potential, etc, and 
have been extensively treated using various field theo- 
retic techniquesi^ii^'^iSii^. The low-energy DOS was 
found to be sensitive to the symmetry of disorders22i^i2^. 
For instance, random vector potential leads the DOS to 
vanish algebraically upon approaching the Fermi surface 
with exponent depending on symmetr y^^i^^ . For this 
kind of disorder, there is essentially no screening effect 
and the Coulomb interaction remains long-ranged, pro- 
vided that the Altshulcr-Aronov type correction to low- 
energy DOS is not included. For random mass potential, 
the zero-energy DOS can have finite value, as a result of 
dynamical discrete symmetry breakin g^^i^^ . In the case 
of weak disorders, the impurity scattering can be treated 
within the conventional self-consistent Born approxima- 
tion, which reveals that the zero-energy DOS acquires a 
finite value of the formic, A^(0) = ;;;^^roln^, with a 

constant scattering rate Fq. The finite A^(0) screens the 
long-range Coulomb interaction. Within the Matsubara 
formalism, such screening effect can be elaborated by in- 
cluding the scattering rate Fq into the polarization func- 
tion. 

To study the role of weak disorders, we first write the 
effective Dirac fermion propagator 




FIG. 3: (a) Dependence of Nc on for different values of T 
at A — > oo; (b) Dependence of Nc on To for different values of 
T at A ^ oo. 



is the same as that appearing in gap equation Eq. ([H]). 
As in the case of chemical potential, the screening effect 
caused by disorder scattering can be directly seen by cal- 
culating the vacuum polarization function 7r(w„,q, Fq) 
and then taking the a;„ — 0, q ^ limit. However, 
even in the instantaneous approximation, it is not easy 
to obtain the complete form of 7r(0, q, Fq). When the 
scattering rate Fq is larger than the thermal scale ^ T, 
Fq > 2ttT, we found that the polarization function can 
be well approximated by the following expression(as de- 
tailed in Appendix) 

N a 

7r(0,q,Fo) « -(q + c'FoCxp(-^)), (15) 
8 c'Fo 



5(iw„,k,Fo) 
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with constant c' = ^i^^^^. At the limit q = 0, it takes 
a finite value 



which contains the scattering rate Fq. Due to the sign 
dependence of scattering rate, the gap equation and po- 
larization function becomes rather complicated. After 
frequency summation within the instantaneous approxi- 
mation, the gap equation is found to be 



m(p,T) 



1 r cP\^ m(k,r) 
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where i){x) is the digamma function. At the clean limit, 
Fq = 0, the imaginary part of digamma function can be 
simplified as Im[-0(i-'--''^^^^£5 



2-kT 
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which is proportional to the scattering rate Fq and de- 
fines the screening factor. Comparing the polarization 
Eq. (|15p with Eq. (7), formally the scattering rate Fq 
plays the role of an effective temperature T. Now the 
interaction function in gap equation Eq. p4p becomes 



mq,ro) 
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q,ro) 
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On the other hand, in the case of small Fq the polariza- 
tion function 7r(0,q, Fq) should be replaced by Eq. (7). 
After solving the full gap equation Eq. (|14p . we present 
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the dependence of Nc on scattering rate Tq for differ- 
ent values of T in Fig. [3]^b) . In order to see the effects 
of screening on gap generation, we also solved the gap 
equation when Tq appears only in the interaction func- 
tion T^(0, q, Tq). The quantitative difference between the 
results in these two cases is negligible. The results in 
Fig. [S^b) show that there is a competition between the 
suppressing effects of thermal fluctuation and disorder 
scattering. At low temperature T, the scattering rate Fq 
dominates; while for small Fq, the thermal effect domi- 
nates. Obviously, a large Fq suppresses the possibility of 
gap generation rapidly. Further, we solved the gap equa- 
tions Eq. ([8]) and Eq. ([9]) with the screening factor simply 
set to be ^ = A^(0) ~ ^j3^Foln^ and found that the 

results are quantitatively similar to Fig. [3Jb). 

One might argue that the low-energy fermionic excita- 
tions are all suppressed once a fermion mass gap opens, 
so the DOS vanishes at energy scale below the gap and 
there is no screening effect. However, for fermion of 
mass m, the zero-energy DOS was found^^ to be N{Q) = 



-Fo In 



In principle, we might include a gap 



m into the polarization function 7r(0,q, Fo,m) and then 
study the gap equation. Since the critical behavior of 
SM-IN transition is studied by linearizing the nonlinear 
gap equation, the mass can be safely set to zero, m — > 0, 
near the bifurcation point. 

Finally, we discuss effect of finite sample volume (area 
in two dimensions). For a graphenc plane of finite spatial 
extent, the particle momenta becomes discrete and the 
momenta transferred in the process of interaction can 
not be arbitrary small. If we still work in the continuum 
field theoretic formalism, this effect can be equivalently 
represented by imposing an infrared cutoff k, given by 
the inverse sample size L~^. Its effects on Nc is nearly 
the same as Fig. [^a) at T = and Fig. ^h) at finite T 
with fi replaced by k, and hence are not shown explicitly. 
The results imply that the sample of large spatial extent 
is more favorable to undergo the SM-IN transitionSi. 

Besides the above four effects, any other mechanism 
that can screen the long-range Coulomb interaction will 
also unavoidably lower the possibility of gap generation. 
If more than one screening effects coexist in reality, the 
suppression of SM-IN transition becomes much more sig- 
nificant, as shown in Fig. [51 In light of these results, 
we conclude that the excitonic insulating state can most 
probably be observed in undoped, clean graphenc of large 
area near absolutely zero temperature. 

Once the long-range Coulomb interaction is screened, 
one interesting question is whether it can be equivalently 
replaced by a short-range or even a contact (on-site) re- 
pulsive interaction^. This question can also be asked in 
another way: is the long-range nature or the strong cou- 
pling nature of Coulomb interaction more important in 
driving the SM-IN transition? If the answer is the lat- 
ter, then the long-range interaction can well be replaced 
by a short-range or contact one. According to the above 
results, it seems that the long-range, rather than strong 



coupling, nature plays the dominant role. As shown in 
Fig. mja), even in the very strong coupling limit A — > c», 
the critical flavor Nc is already less than the physical 
flavor 2 when the screening factor fi ~ 10~^. For mod- 
erately strong coupling A = 2.5, the excitonic insulating 
behavior becomes impossible even if the screening factor 
is only as small as fi ^ 10" 
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FIG. 4: Dependence of Nc on fi for different g. 

In order to test the role of contact interaction and see 
its difference from the screened Coulomb interaction, we 
add one quartic interacting term to the Hamiltonian. 
There are several choices for the four-fermion coupling 
term, classified by the gamma matrices used to define 
the actio n^^i^^ . For simplicity, we consider only one of 
them, i.e., 
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N 
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(18) 



To the lowest order, this contact interaction contributes 
the following term 
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rf^k m(k, /?) tanh 



2T 



A/k2 + m2(k,/3) 



(19) 



to the gap equation, where the dimensionless coupling 
is g = NGA/vp and the scaling Wi?k — > k is made as 
before. The whole gap equation is solved with results 
shown in Fig. Hat T = 10-^cV{^ lOmK). The contact 
four-fermion interaction has opposite effect on the crit- 
ical fiavor Nc as compared with the screening factor /z: 
while the latter rapidly suppresses Nc, the former is very 
efficient in promoting the system towards the excitonic 
insulating phase (note there is no Goldstone boson in the 
insulating phase since the total Hamiltonian preserves 
discrete chiral symmetry ip 75 V')- Thus we see that the 
contact four-fermion interaction is actually different from 
the screened Coulomb interaction. For a relatively large 
screening factor ^, the latter is unable to generate exci- 
tonic gap even in the A — *■ 00 limit, while the former can 
generate such gap when its coupling is larger than some 
critical value g > gc- The reason for this difference can 
be seen from the gap equation: for screened Coulomb in- 
teraction, q appearing in the denominator suppresses the 
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contribution from large momenta, while fi in the denom- 
inator suppresses the contribution from small momenta; 
on the contrary, for contact fermion interaction, the cou- 
pling g is constant in the whole momenta region without 
any suppressing effect. In conclusion, the SM-IN transi- 
tion is still possible if there is additional strong contact 
fermion interaction, even when the screened Coulomb in- 
teraction itself can not open the gap. 

We end with a brief discussion on the validity of the 
gap equation used in this paper. In a rigorous treatment, 
the excitonic gap generation should be studied by solv- 
ing the self-consistent equations of fermion self-energy 
function, wave function renormalization. Coulomb inter- 
action, and vertex function. In practice, a number of 
approximations must be utilized. Here we kept only the 
Fock diagram for the fermion self-energy and omit all 
higher order corrections of the expansion^. The 
results should be qualitatively reliable for large N. To 
verify the conclusions obtained in the leading order, it 
would be necessary to include these corrections (such 
as wave function renormalization, vertex function cor- 
rection, etc.) since they might change the quantities of 
critical parameters Ac and fic considerably for the physi- 
cal flavor N = 2. However, this is beyond the scope of the 
present work. Another question concerns the important 
effect of velocity renormalization on the excitonic gap 
instabilitji^ii^. This effect has been addressed recently 
by incorporating the momentum-dependent fermion ve- 
locity into the gap equation^l. It was found that the 
velocity renormalization does not dramatically affect the 
excitonic instability^!. 
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APPENDIX A: CALCULATION OF 7r(0, q, /?, To) 

In this appendix, we present the derivation of the po- 
larization function at finite impurity scattering rate and 
finite temperature. The fermion contribution to the vac- 
uum polarization is given by 



d^k Tr[7o^7o(^ + ^)] 



)2 k^{q + ky 



(Al) 



Here qq = iujm = and fco = iuJn = (2n+i)£^ ^ 
new momentum variable is defined hy I = k + xq with 

Within the instantaneous approximation = 0, the 
polarization function reduces to 

7r(0, q, /?) - ^ dx J ^ [Si - 21^82] , (A2) 
where 5^=1^2 is given by 

5, = y (A3) 

^t-^iii+p+xii-xwr 

In the presence of impurity scattering rate Fg, the vari- 
able lo should be replaced by 

l'0^^{n+l + ^FoSgnw„). (A4) 



Using th e notation in Re fi^, we define a new variable 
Y = ^ + x{l — x)q2. Using the identity 



SiX,Y) = J2 



{n + a:)2 + y2 



— [^{X + tY)-^{X-tY)], (A5) 



the function Si now becomes 



^ + ^ro,y) + 5(i-i + Aro,y) 



= ^lmml+X' + ^Y)], (A6) 

from which the function 6*2 is given by ^2 = 
Sy ^- Define t = fY, and ^ = + C^] = 
-f (■y/a;(l — a;)q2)2]. then the polarization function 
7r(0, q, /?, Fq) is written as the following integral 



, , , 

— — I dx dt 



.-'0 



t^-C^ d , ,,1 13 /3 



(A7) 



It is hard to compute this integral analytically. For rela- 
tively large scattering rate Fg, we found that the "0 func- 
tion can be approximated by the analytic expression 



,1/3 /3 , TT nt 

ip(- + —To + i—t) « - tanh — r-^ — 
^^2 2tt ' 2n ' 2 lO'^^p 



(A8) 



for > 1 (with error 1% for 27 To ^ 1 and averag- 

ing error 5% for ^Fg « 1). Then the integration over 
variable t can be carried out with the result 



^(0,q,Fo)«^^ dx-ln[2cosh ^ ]■ 



(A9) 



It has the similar form as Eq. ^ with Fg playing the 
role of an effective "temperature", thus the polarization 
function can now be approximated by 

N a 
7r(0, q, Fg) « ^(q + c'Fg cxp(-^)), (AlO) 

o C i g 

where c' = ^'"^2°'"^ ■ At the clean hmit, Fg < n/P, we 
have 

In this case, the polarization function is still approxi- 
mated by Eq. ([7]). 
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